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Various modern nucleon-nucleon (NN) potentials yield a 
very accurate fit to the nucleon-nucleon scattering phase 
shifts. The differences between these interactions in describ- 
ing properties of nuclear matter are investigated. Various 
contributions to the total energy are evaluated employing the 
Hellmann - Feynman theorem. Special attention is paid to the 
two-nucleon correlation functions derived from these interac- 
tions. Differences in the predictions of the various interactions 
can be traced back to the inclusion of non-local terms. 



I. INTRODUCTION 

The microscopic theory of nuclear structure based on 
realistic nucleon-nucleon (NN) interactions is a very de- 
manding subject because it requires the description of a 
strongly correlated many-fermion system. Attempts to 
determine e.g. the energy of nuclei from a realistic NN 
interaction by using the mean-field or Hartree-Fock ap- 
proximation fail badly: such attempts typically yield un- 
bound nuclei. The strong short-range and tensor compo- 
nents of a realistic NN interaction induce correlations into 
the many-body wavefunction of nuclear systems. Many 
attempts have been made to measure these correlations 
in detail. As an example for such measurements we men- 
tion the exclusive (e, e'NN) reactions, which have been 
made possible using modern electron accelerators [Q,^. 
The hope is that the detailed analysis of such experi- 
ments yields information about the correlated wavefunc- 
tion of the nucleon pair absorbing the virtual photon. 
This could be a very valuable test for the model of the 
NN interaction producing these correlations. 

In recent years, a new generation of realistic NN po- 
tentials has been developed which produce very accurate 
fits of the proton-proton and proton-neutron (pn) scat- 
tering phase shifts |^-^ . Since these fits are based on the 
same phase shift analysis by the Nijmegen group Q and 
yield a value for the x^/datum very close to one, these 
various potentials could be called phase-shift equivalent 
NN interactions. 

This means that the on-shell matrix elements for the 
transition matrix T are essentially identical. This, how- 
ever, does not imply that the underlying potentials nor 
the effective interaction between off-shell nucleons mov- 
ing inside a nucleus are the same. Indeed it has been 



demonstrated that these phase shift equivalent potentials 
yield different results even for the deuteron. Of course 
all of them reproduce the same empirical binding energy 
and other observables, because these are part of the ob- 
servables to which the interaction has been fitted. How- 
ever, the various contributions to the total energy, the 
kinetic energy and the potential energy in the ^Si and 
'^Di partial waves, are quite different indicating that also 
the two-body wavefunctions must be different j7|J|] . 

It is one of the aims of this study to explore whether 
similar differences can also be observed in calculating the 
energy of nuclear matter. These energies are calculated 
using the Brueckner-Hartree-Fock (BHF) approximation. 
The BHF approach for nuclear matter assumes a model- 
wavefunction of a free Fermi gas, occupying plane wave 
states up to the Fermi momentum. The effects of corre- 
lations are taken into account in terms of the Brueckner 
G-matrix. The BHF approximation does not give direct 
access to quantities like the kinetic or the potential en- 
ergy. However, in the next section we will illustrate how 
the Hellman-Feynman Theorem [|| can be used to calcu- 
late these quantities and also the expectation value of the 
TT exchange calculated for the correlated wave function. 
The results for various contributions to the binding en- 
ergy as well as the wave functions of correlated NN pairs 
in nuclear matter will be presented in section 3. The last 
section contains a summary and conclusions. 



II. CORRELATIONS IN BRUECKNER HARTREE 
FOCK 

The central equation of the BHF approximation is the 
Bethe-Goldstone equation, which defines an effective in- 
teraction G for two nucleons in nuclear matter occupying 
the plane wave states i and j by 



G\ij >= V\ij > +V- 
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Ho 



(1) 



Here and in the following V stands for the bare NN in- 
teraction, Q denotes the Pauli operator, which prevents 
of the interacting nucleons into intermediate states with 
momenta below the Fermi momentum kp, Hq defines 
the spectrum of intermediate two-particle states and the 
BHF single-particle energies are defined by 
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as the sum of the kinetic energy of a free nucleon with 
mass TO and momentum ki and the potential energy. The 
single particle potential corresponds to the Hartree-Fock 
approximation but calculated in terms of the effective 
interaction G rather than the bare interaction V. Also 
the total energy of the system is calculated in a similar 
way containing the kinetic energy per nucleon of a free 
Fermi gas 
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and the potential energy calculated in the Hartree-Fock 
approximation replacing V by the effective interaction 
G (For a more detailed description see e.g. 0). This 
means that the BHF approach considers a model wave 
function, which is just the uncorrelated wave function of 
a free Fermi gas and all information about correlations 
are hidden in the effective interaction G. Since this ef- 
fective interaction is constructed such that G applied to 
the uncorrelated two-body wave function yields the same 
result as the bare interaction V acting on the correlated 
wave function 



G\ij >= V\ij >corr. 



(4) 



the comparison of this equation with (||) allows the defi- 
nition of the correlated two-nucleon wave function as 



.= I«J > 
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This representation demonstrates that the correlated 
wave function contains the uncorrelated one plus the so- 
called defect function, which in this approach should drop 
to zero for relative distances between the two nucleons, 
which are larger than the healing distance. 

The BHF approach yields the total energy of the sys- 
tem including effects of correlations. Since, however, it 
does not provide the correlated many-body wave func- 
tion, one does not obtain any information about e.g. the 
expectation value for the kinetic energy using this corre- 
lated many-body state. To obtain such information one 
can use the Hellmann-Feynman theorem, which may be 
formulated as follows: Assume that one splits the total 
Hamiltonian into 



H = Ho + AV 



(6) 



and defines a Hamiltonian depending on a parameter A 
by 



(A) =Ho + XAV . 
If E\ defines the eigenvalue of 



(7) 
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the expectation value of AV calculated for the eigen- 
states of the original Hamiltonian H = H{1) is given as 



(9) 
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The BHF approximation can be used to evaluate the 
energies Ex, which also leads to the expectation value 
< \E'|AT^|^' > employing this eq.(||). In the present work 
we are going to apply the Hellmann-Feynman theorem 
to determine the expectation value of the kinetic energy 
and of the one-pion-exchange term AV = Vtt contained 
in the different interactions. 



III. RESULTS AND DISCUSSION 

The main aim of the work presented here is to in- 
vestigate differences in nuclear structure calculations 
originating from four different realistic NN interactions, 
which are phase-shift equivalent. These four interactions 
are the so-called charge-dependent Bonn potential (CD- 
Bonn) [pj , the Argonne V18 (ArV18) D and the versions 
I (Nijmlj and H (Nijm2) of the Nijmegen interaction 
All these models for the NN interaction include a one- 
pion exchange (OPE) term, using essentially the same 
ttNN coupling constant, and account for the difference 
between the masses of the charged (7r±) and neutral (ttq) 
pion. However, even this long range part of the NN 
interaction, which is believed to be well understood, is 
treated quite differently in these models. The Nijmegen 
and the Argonne V18 potentials use the local approxi- 
mation, while the pion contribution to the CDBonn po- 
tential is derived in a relativistic framework assuming 
pseudoscalar coupling. It has recently been shown that 
the non-localities included in the relativistic description 
of the CDBonn potential tends to lead to smaller D-state 
probabilities in the deuteron ^ . 

The description of the short-range part is also differ- 
ent in these models. The NN potential Nijm2 is a 
purely local potential in the sense that it uses the local 
form of the OPE potential for the long-range part and 
parameterizes the contributions of medium and short- 
range distances in terms of local functions (depending 
only on the relative displacement between the two inter- 
acting nucleons) multiplied by a set of spin-isospin op- 
erators. The same is true for the Argonne Vis potential 
[Q. The NN potential denoted by Nijml Q uses also 
the local form of OPE but includes a term in the 
medium- and short-range central-force (see Eq. (13) of 
Ref. [^) which may be interpreted as a non-local con- 
tribution to the central force. The CD-Bonn is derived 
in the framework of the relativistic meson field theory. 
It is calculated in momentum space and contains non- 
local terms in the short-range as well as long-range part 
including the pion-exchange contribution. 

First differences in the prediction of nuclear properties 
obtained from these interactions are displayed in table || 
which contains various expectation values calculated for 
nuclear matter at the empirical saturation density, which 
corresponds to a Fermi momentum kp of 1.36 fm^^. The 
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most striking indication for the importance of nuclear 
correlations beyond the mean field approximation may 
be obtained from the comparison of the energy per nu- 
cleon calculated in the mean-field or Hartree-Fock (HF) 
approximation. All energies per nucleon calculated in 
the (HF) approximation are positive, therefore far away 
from the empirical value of -16 MeV. Only after inclusion 
of NN correlations in the BHF approximation results are 
obtained which are close to the experiment. While the 
HF energies range from 4.6 MeV in the case of CDBonn 
to 36.9 MeV for Nijm2, rather similar results are obtained 
in the BHF approximations. This demonstrates that the 
effect of correlations is quite different for the different in- 
teractions considered. However it is worth noting that all 
these modern interactions arc much "softer" than e.g. the 
old Reid soft-core potential in the sense that the HF 
result obtained for the Reid potential (176 MeV) is much 
more repulsive. 

Another measure for the correlations is the enhance- 
ment of the kinetic energy calculated for the correlated 
wave function as compared to the mean field result which 
is identical to Tpc the energy per particle of the free 
Fermi gas. At the empirical density this value for Tpc is 
23 MeV per nucleon. One finds that correlations yield an 
enhancement for this by a factor which ranges from 1.57 
in the case of CDBonn to 2.09 for Nijml. It is remark- 
able that the effects of correlations, measured in terms of 
the enhancement of the kinetic energy or looking at the 
difference between the HF and BHF energies, are signif- 
icantly smaller for the interactions CDBonn and Nijml, 
which contain non-local terms. 

The table | also lists the expectation value for the pion- 
exchange contribution Vt^ to the two-body interaction. 
Here one should note that the expectation value of V,r 
calculated in the HF approximation is about 15 MeV al- 
most independent of the interaction considered. So it 
is repulsive and completely due to the Fock exchange 
term. If, however, the expectation value for 14- is evalu- 
ated for the correlated wave function, one obtains rather 
attractive contributions ranging from -22.30 MeV per nu- 
cleon (CDBonn) to -40.35 MeV (ArV18). This expecta- 
tion value is correlated to the strength of the tensor force 
or the D-state probability Pd calculated for the deuteron 
(see table ||as well). Interactions with larger Po, like the 
ArVlS, yield larger values for < 1/^ >. For a further sup- 
port of this argument we also give the results for three 
different version of charge-independent Bonn potentials 
A, B and C, defined in |ll|| . 

All this demonstrates that pionic and tensor correla- 
tions are very important to describe the binding proper- 
ties of nuclei. In fact, the gain in binding energy due to 
correlations from T4- alone is almost sufficient to explain 
the difference between the HF and BHF energies. 

Until now we have just discussed results for nuclear 
matter at one density. The values for the kinetic energy, 
< T >, and < > are displayed for various densities 
in Fig. 0. One finds that the ratio of the kinetic energy 
calculated for the correlated wave function, < T >, and 



the energy of the free Fermi-gas < Tpc > decreases as 
a function of density. This plot furthermore shows that 
the results for the different interactions can be separated 
in two groups: the local interactions, ArV18 and Nijm2, 
yield larger kinetic energies than CDBonn and Nijml, 
which contain nonlocal terms. 

The lower part of Fig. |l| shows that the pionic con- 
tribution to the total energy is quite different for the 
interactions. It is strongest for ArV18, getting more at- 
tractive for larger densities. The pionic contribution ob- 
tained from the other potentials is weaker and does not 
exhibit this increase at high densities. This may indicate 
that the enhancement of pionic correlations, which has 
been discussed in the literature as an indication for pion 
condensation ||l3|] , is a feature which may not be repro- 
duced by realistic interactions different from the Argonne 
potentials. 

A different point of view on nuclear correlations may be 
obtained from inspecting the the relative wave functions 
for a correlated pair \ij >corr. defined in (||). Results 
for such correlated wave functions for a pair of nucle- 
ons in nuclear matter at empirical saturation density are 
displayed in Figs |^ and ^j. As an example we consider 
wave functions which "heal" at larger relative distances 
to an uncorrelated two-nucleon wave function with mo- 
mentum q — 0.96 fm^^ calculated at a corresponding 
average value for the starting energy. 

Fig. 1^ shows relative wave functions for the partial 
wave ^5*0. One observes the typical features: a reduc- 
tion of the amplitude as compared to the uncorrelated 
wave function for relative smaller than 0.5 fm, reflecting 
the repulsive core of the NN interaction, an enhancement 
for distances between « 0.7 fm and 1.7 fm, which is due 
to the attractive components at medium range, and the 
healing to the uncorrelated wave function at large r. One 
finds that the reduction at short short distances is much 
weaker for the interactions CDBonn and Nijml than for 
the other two. This is in agreement with the discussion 
of the kinetic energies (see Fig. |]) and the difference be- 
tween HF and BHF energies (see table |). The nonlocal 
interactions CDBonn and Nijml are able to fit the NN 
scattering phase shifts with a softer central core than the 
local interactions. 

Very similar features are also observed in the ^Si par- 
tial wave displayed in the left half of Fig. For the '^Di 
partial wave, shown in the right part of Fig. ^, one ob- 
serves a different behavior: All NN interactions yield an 
enhancement of the correlated wave function at r « 1 
fm. This enhancement is due to the tensor correlations, 
which couples the partial waves "^^i and '^Di. This en- 
hancement is stronger for the interactions ArVlS, Nijml 
and Nijm2 than for the CDBonn potential. Note that 
the former potential contain a pure nonrelativistic, local 
one-pion-exchange term, while the CDBonn contains a 
relativistic, nonlocal pion-exchange contribution. 

This behavior in the coupled '^Si and ^Di waves can 
also be observed in the corresponding wave functions for 
the deuteron, plotted in Fig. ^. 
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IV. CONCLUSIONS 

Four modern NN interactions, the charge-dependent 
Bonn potential (CDBonn), the Argonne V18 (ArV18) 
and two versions of the Nijmegen potential (Nijml and 
Nijni2), which all give an excellent fit to NN scattering 
phase shifts, exhibit significant differences in calculating 
NN correlation functions and other observables in nuclear 
matter. Two of these interactions, CDBonn and Nijml, 
contain nonlocal terms. These two interactions are con- 
siderably softer than the other interactions. This con- 
clusion can be derived from three different observations: 
The Hartree-Fock energies are less repulsive, the kinetic 
energies calculated with the correlated wave functions 
are smaller and the correlated wave function in relative 
S states are less suppressed at small relative distances. 

The interactions also differ quite significantly in the 
pionic or tensor correlations they induce. This is indi- 
cated to some extent by the deuteron wave function, in 
particular by the D-state probability. These differences, 
however, are even enhanced in the nuclear wave functions 
leading to drastic differences in the pionic contribution 
to the nuclear binding energy. The Argonne potential 
in particular yields a large pionic contribution, which in- 
creases with density. This importance of the pionic cor- 
relations is not observed for the other interactions. 

It would be of great interest to study whether the dif- 
ferences between the correlations predicted from these 
interactions can be observed in experiments like the ex- 
clusive (e, e'iViV) reactions in order to discriminate the 
various models for the NN interaction. 
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CDBonn 


ArV18 


Nijml 


Nijm2 


A 


B 


C 


Reid 


<E> 


-17.11 


-15.85 


-15.82 


-13.93 


-16.32 


-15.32 


-14.40 


-12.47 


<V > 


-53.34 


-62.92 


-55.08 


-61.94 


-52.44 


-53.03 


-54.95 


-61.51 


<T> 


36.23 


47.07 


39.26 


48.01 


36.12 


37.71 


40.55 


49.04 


< K > 


-22.30 


-40.35 


-28.98 


-28.97 


-12.48 


-26.87 


-45.74 


-27.37 


<i5>HF 


4.64 


30.34 


12.08 


36.871 


7.02 


10.07 


29.56 


176.25 


Pd [% ] 


4.83 


5.78 


5.66 


5.64 


4.38 


4.99 


5.62 


6.47 



TABLE I. Energies calculated for nuclear matter with Formi momentum kp = 1.36 fm^^. Results are listed for the energy 
per nucleon calculated in BHF (< E >) and Hartree-Fock (< E >hf) approximation. Furthermore the expectation value for 
the NN interaction < V >, the kinetic energy < Tj^jjj > and the one-pion-exchange term < 14- > are listed. For completeness 
we also give the D-state probability calculated for the deuteron Pn ■ Results are presented for the charge-dependent Bonn 
(CDBonn) the Argonne V18 (ArV18) Q and two Nijmegen (Nijml, Nijm2) Q interactions. For a comparison results 
are also given for three older versions of the Bonn interaction (A,B,C) |l^ and the Reid soft core potential ||l^, which is 
supplemented in partial waves in which it is not defined by the OBE C potential. All energies are given in MeV per nucleon. 
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FIG. 1. The upper part of this figure displays the ratio of the kinetic energy per nucleon, calculated for the correlated 
state, to the energy per nucleon of a free Fermi-gas ^ as a function of the Fermi momentum. The lower part exhibits the 
expectation value of the one-pion-exchange contribution to the binding energy per nucleon. Different realistic NN interactions 
are considered. 
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FIG. 2. Correlated wave functions \ij >corr. as defined in as a function of the relative distance for the ^So partial 
wave. Results are shown for a pair of nucleons in nuclear matter at empirical saturation density, which heal to an uncorrelated 
two-nucleon wave function with momentum q — 0.96 fm— 1 at larger distances. The curves are labeled by the interactions, 
which were considered. 
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FIG. 3. Correlated wave functions as a function of the relative distance for the ^Si and ^Di partial waves. Further details 
see Fig. |. 
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